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ABSTRACT 

It has been observed that hot Jupiters located within 0.08 AU of their host stars commonly display 
radii in excess of those expected based on models. A number of theoretical explanations for this 
phenomenon have been suggested, but the ability of any one mechanism to account for the full range 
of observations remains to be rigorously proven. I identify an additional heating mechanism, arising 
from the interaction of the interplanetary magnetic field and the planetary magnetosphere, and show 
that this is capable of providing enough energy to explain the observed planetary radii. Such a 
model predicts that the degree of heating should be dependent on the stellar magnetic field, for which 
stellar activity serves as a proxy. Accordingly, I examine populations of hot Jupiters from the Kepler 
database and confirm that stellar activity (determined using Kepler CDPP levels) is correlated with 
the presence of planetary radii inflated beyond the basal level of i? = 0.87Rj identified by previous 
researchers. I propose that the primary mechanism for transferring energy from the magnetosphere 
to the planetary interior is Joule heating arising from global electric circuits analogous to those seen 
in solar system objects. 

Subject headings: stars: planetary systems — planets and satellites: general — planets and satellites: 
atmospheres 



1. INTRODUCTION 

Giant planets are intrinsically interesting both because 
they enable us to better understand the formation and 
evolution of our own solar system and because they serve 
as laboratories for the science of materials at tempera- 
tures and pressures beyond those achievable in the labo- 
ratory (Perna, Heng, & Pont 2012). The study of extra- 
solar giant planets became an observational discipline in 
1995, with the discovery of the first example, 51 Peg b 
(Mayor & Queloz 1995). Four years later, Charbonneau 
et al. (2000) found the first transiting extrasolar giant 
planet, HD 209458b. The existence of transits enabled 
measurement of the planetary radius, which Charbon- 
neau et al. found to be inflated beyond both expectations 
based on models and those based on our experiences in 
our own solar system. Since that time, approximately 
200 hot Jupiters have been discovered, and inflated radii 
have proven to be a common, although not universal, 
phenomenon (Demory & Seager 2011). 

Planets in this size range are supported by a combi- 
nation of electron degeneracy and thermal pressure, so 
the resulting models predict an approximately constant 
radius for planets with masses between ~ 1 and ~ 7Mj. 
Quite early on, Guillot et al. (1996) noted that, while gi- 
ant planets far from their host stars tend to follow the ex- 
pected mass-radius relation, those closer in are typically 
inflated in comparison to model predictions. The most 
recent study in this area is by Demory & Seager (2011), 
who examined a Kepler sample of 138 giant planet can- 
didates, and determined that those experiencing incident 
stellar fluxes less than ~ 2 x 10^ erg cm~^ s~^ showed 
constant radii of 0.87 ± 0.12Rj, while planets experienc- 
ing stellar fluxes above this level were inflated, with the 
degree of inflation increasing along with the stellar flux 
level. 
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Several models have been put forward to explain these 
observations. A more detailed explication is given in 
Fortney & Nettlemann (2010), but essentially most mod- 
els require either modiflcations to improve heat retention, 
or additional energy sources. Thus, Burrows et al. (2007) 
suggest enhanced atmospheric opacities, while Chabrier 
& Baraffe (2007) show that heavy element abundance 
gradients could suppress convection, both leading to 
greater radii. Alternatively, Hansen & Barman (2007) 
suggest that preferential He mass loss would lead to H- 
rich compositions, and thus greater radii. In the second 
category, additional energy sources have been posited 
to arise from tidal heating (Bodenheimer et al. 2001; 
Gu, Bodenheimer, & Lin 2004), enhanced conversion of 
the incident stellar flux to planetary wind energy (Guil- 
lot & Showman 2002), or interaction of ionized plane- 
tary winds with the planetary magnetic held (Batygin & 
Stevenson 2010). The amount of additional heating re- 
quired can be as large as 8 x 10^^ erg s~^ (for Tres-4b; see 
Batygin & Stevenson 2010), and the additional heating 
must be supplied to the deep interior of the planet, typ- 
ically below the 10^ — 10'^ bar level. Furthermore, some 
proposed models have difficulties in accounting for the 
increase in radii for large incident stellar fiuxes (Laugh- 
lin et al. 2011), although Batygin et al. (201 1) show 
that the "Ohmic dissipation" model, whereby heating 
derives from interactions between atmospheric flows and 
the planetary magnetic field, can successfully reproduce 
this phenomenon. 

The critical incident stellar flux level of ~ 2 x 
10* erg cm~^ s~^ identified in Demory & Seager trans- 
lates, in the case of a star like the Sun, to a circular 
orbit with a semimajor axis of approximately 0.08 AU, 
or 18Ro. It is suggestive that this roughly corresponds 
to the Alfven radius for the Sun, variously estimated at 
5 - 5OR0 (Aibeo et al. 2007, Scrijver, DeRosa, & Title 
2003). In our solar system, magnetospheric extraplane- 
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tary energy input at both the Earth and Jupiter is dom- 
inated by kinetic energy from the solar wind, but in this 
paper I suggest that for the much closer hot Jupiters the 
magnetic energy component of the stellar wind is a more 
important contributor and may be capable of providing 
enough energy to account for the inflated radii which are 
observed for hot Jupiters close to their host stars. 

2. MODEL 

Any planet with an internal magnetic field can be 
expected to be surrounded by a magnetosphere, which 
serves as an energy reservoir, with energy flowing in from 
the stellar wind and the planetary rotation, and out via 
heating of the solar wind and planet through a number 
of mechanisms (Vasyliunas 2010). In this work, I ne- 
glect the influence of planetary rotation, and consider a 
planet with radius RpL ~ Rj and a magnetic dipole field 
of 1 - 25G. 

Pressure balance considerations then yield the scale 
size of the magnetosphere. Stellar wind pressure has 
contributions from the stellar magnetic pressure Pb = 
i?^/87r, the kinetic (ram) pressure Pk = /0Wu,^/2, and 
thermal pressure Pth = nkT] the latter can be neglected 
as it is at least two orders of magnitude less than the 
other contributions over the entire range of interest. In 
the case of both the Earth and Jupiter, the ram pressure 
term dominates over stellar magnetic pressure, but the 
presence of multipole components in the stellar magnetic 
field implies that this may not be the case for small or- 
bital radii. Balancing stellar wind pressure against the 
pressure from the planetary dipole magnetic field gives 

Rm = RpL { 2 , /" — 2 r 

This is a generalization of the Chapman-Ferraro radius 
(Chapman & Ferraro 1931), such that the pressure of the 

stellar wind is applied over the area a = ttRm"^ ■ 

I model the stellar wind using the model outlined in 
Suzuki (2006), with the magnetic field given by the ana- 
lytic description in Banaszkiewicz, Axford, & McKenzie 
(1998). This relatively simple model gives results similar 
to those seen in more complex implementations for radii 
outside the source surface (taken to be Rss = 2.5R,), 
and in any case results here are insensitive to details of 
the model used. The calculated field {Bp, Bz) is pre- 
sumed to possess cylindrical symmetry; in addition, I 
limit myself to the magnetic field and wind in the equa- 
torial plane of the star (presumed to be the orbital plane 
of the planet), but the model is easily extended to more 
complex (and realistic) cases. I followed Banaszkiewicz, 
Axford, & McKenzie in adopting if = 1.0, M = 1.789, 
ai — 1.538, and Q = 1.5 in order to force the wind model 
into close correspondence with the observed solar wind. 

I adopt a commonly-used empirical formula describing 
the rate of terrestrial magnetospheric heating by storm 
and sub-storm processes (the so-called "epsilon parame- 
ter", see, e.g. Stern 1984): 

e = vIq^B'^ sin^ 6 (2) 

Here 9 is the polar angle of the solar interplanetary mag- 
netic field vector, projected onto the Y-Z plane (Akasofu 
1981). Such a model can be defended on dimensional 
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Figure 1. Contributions to tiie total stellar wind energy available 
to the planetary magnetosphere. The kinetic energy component is 
shown in red (dot-dash), while the magnetic contribution is in blue 
(dotted); total wind energy is given in black (solid). The thermal 
energy component is not shown as it is more than two orders of 
magnitude smaller than the others over the range of the figure. 
Note that for semimajor axes below about 30Rq , corresponding to 
approximately 0.15 AU, the magnetic contribution dominates. In 
all figures, the abscissa runs from a = Rss to a = 1 AU. 

grounds, as well as in more detailed arguments (Vasyli- 
unas 2010), which identify the parameter as the power 
generated by the dynamo defined by the interaction of 
solar wind and magnetosphere. Note that Iq is a scale 
length for the magnetosphere. Early work on the ter- 
restrial magnetosphere adopted Iq = 1R@, but more re- 
cent analysis (Tanskanen et al. 2005) has adopted a re- 
vised scaling of Iq = 10i?e- Noting that (unsurprisingly) 
Iq ~ Rm for both the Earth and Jupiter, in this work I 
adopt Iq = Rm. 

3. RESULTS 

Figure 1 illustrates the rising importance of the mag- 
netic energy component of the stellar wind pressure for 
smaller semi-major axes. For orbital distances below 
about 0.15 AU, it is dominant. Note in particular that 
at a = 0.08AU, the available wind magnetic power is ap- 
proaching 10^^ erg s~^ and rising rapidly; this is just the 
point at which Demory & Seager (2011) see hot Jupiter 
radii start to increase. In addition, the power supplied by 
the wind is of the correct order of magnitude to supply 
the "missing" planetary heating component. 

Of course, this agreement could simply be a fortu- 
itous consequence of well-chosen values for the two main 
free parameters, -BpL and B*. Decreasing the planetary 
magnetic dipole strength BpL leads to smaller planetary 
magnetospheres, and thus to less magnetospheric energy 

deposition. However, since as a — > RsSi Rm oc Bp^^, the 
effect is quite small, and thus the energy deposition seen 
in the model is relatively insensitive to Bp^. 

The effects of changing the stellar magnetic field are 
much more substantial, as shown in Figure 2, which illus- 
trates the effects of increasing the stellar magnetic field 
by factors of 10 and 100, while maintaining all other pa- 
rameters constant. Cranmcr & Saar (2011) suggest that 
the photospheric magnetic filling factor is related to the 
stellar Rossby number by /* ~ Ro~", where n lies be- 
tween 2.5 and 3.4. If the magnetic field at Rss scales 
similarly, this would imply rotation rates of 4 — 6 days 
and filling factors near saturation for the most extreme 
example shown. The magnetic component of wind power 
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Figure 2. The effects of changing the stellar magnetic field. The 
solid (black) line shows the solar case, as in Figure 1, while the dot- 
dash (red) and dotted (blue) lines show the impact of increasing 
the stellar magnetic field by factors of 10 and 100, respectively. 
Even in the solar case, 10^ erg is available at a = 0.1 AU, 
while in the most active case, energy deposition rates as high as 
IqSO gj.g g— 1 g^j.g possible at a = 0.05 AU. 

which is available to the magnetosphere in this case can 
be as large as 10^° erg s~^ at a = 0.04AU. 

In cases where the stellar magnetic field is substan- 
tially greater than the solar case, results are essentially 
independent of mass loss rate, for any mass loss rate 
which is reasonably solar-like. This is because, for such 
stars, magnetic effects are dominant in determining both 
Rm and wind power available to the planetary magneto- 
sphere for orbits inside 0.08 AU. 

4. DISCUSSION 

4.1. Observational Testing 

The magnetic heating model proposed above has ob- 
servable ramifications, which can be tested using the rel- 
atively large and self-consistent data set available from 
the Kepler mission. In particular, as seen in Figure 2, the 
magnetospheric heating rate in this model is strongly de- 
pendent on stellar magnetic field strength, encapsulated 
at photospheric levels in the parameter -B*/, . In turn, 
i?*/* is closely related to stellar activity levels (Pizzolato 
et al. 2003), for which typical observational proxies are 
normalized soft X-ray luminosity or emission in the cores 
of the Ca II H and K lines (Hah 2008). Ahhough these 
particular proxies are generally unavailable for most stars 
with hot Jupiters (only 11 of the hot Jupiter host stars 
have extant Ca II measurements), I will follow Basri et al. 
(2011) and others in adopting the Kepler Combined Dif- 
ferential Photometry Precision (CDPP) as a measure of 
stellar photometric variability, and hence activity. A first 
order attempt was made to remove photometric noise 
from the CDPP by subtracting (in an rms sense) the 
simple model CDPP derived in Machalek et al. (2011, 
see also Gilliland et al. 2011). Corrections were minor 
in most cases, and do not substantially affect my conclu- 
sions. 

Approximately following Demory & Seager, I 
used the Kepler Planet Candidate Data Explorer 
([http: / / planetquest.jpl.nasa.gov /kepler) to select all 
KUl planet candidates with HK^ < i?PL < 22i?0, 
discarding those with no reported value for CDPP, 
resulting in 153 candidates. Figure 3 shows the distri- 
bution of CDPP as a function of planetary radius, with 
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Figure 3. Corrected CDPP values as a function of planetary ra- 
dius for hot Jupiters (8i?0 < RpL < 22/?0) selected from KOI. 
CDPP values, corrected for non-astronomical noise sources, serve 
as a proxy for stellar activity. The vertical red line corresponds to 
the constant radius of 0.87 it 0.12Rj found by Demory & Seager 
(2011) for low incident stellar fiux, and the dotted lines show the 
ztlcr levels. Visually, the populations vary at about that point, in 
the sense that while large planetary radii are associated with both 
active and inactive stars, small planetary radii are only associated 
with inactive ones. 

the constant radius level of 0.87Rj seen by Demory & 
Seager (2011) indicated by the red vertical line. While 
large planetary radii are associated with both active and 
inactive stars, small planetary radii are apparently only 
associated with inactive stars. I quantified this observa- 
tion by applying a Kolmogorov-Smirnov test (Hollander 
& Wolfe 1973) to test the hypothesis that the two groups 
represent the same underlying population. The K-S test 
indicates a probability p = 0.03 that the hypothesis is 
true, so we can conclude that a correlation exists be- 
tween inflated radii and stellar activity. However, there 
is clearly considerable scatter in such a correlation (see 
Figure 3), as one might expect considering the impact 
of other parameters such as planetary composition, the 
strength and morphology of the planetary magnetic 
field, noncircular orbits, varying orbital alignments, 
and possible variations in the efficiency of wind energy 
deposition into the planetary interior. 

4.2. Possible Mechanisms 

We expect that, as in our solar system, exoplanet mag- 
netospheres will be electrically connected to the upper 
ionosphere via a system of field aligned currents (FAC), 
the net result of which is to map the open solar wind mag- 
netic field lines to the polar ionosphere and produce an 
electric potential structure across the ionosphere (Lysak 
1980, Zhang et al. 2011). The scale of the potential dif- 
ference can be estimated from AV 2r]VwBwRM^ where 
ry is a scale factor of order unity. In the solar system, tak- 
ing 77 = 1 gives Ay® ~ 180 kV and AVj ~ 2.2 MV, both 
of which are of the correct order of magnitude compared 
to observations (Singh et al. 2007, Cowley et al. 2003). 
For a hot Jupiter located at a = 0.08 AU, the resulting 
potential difference is AV = 320 MV. 

From the ionosphere, energy can be transferred down- 
wards by a number of mechanisms (Vasyliunas, 2010), 
but here I will consider only two. First I consider auroral 
particle precipitation, in which relativistic particles are 
accelerated downwards from the ionosphere. While the 
details of the acceleration mechanism are complex and 
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somewhat obscure (Olsson & Janhunen, 2003), an order- 
of-magnitude energy estimate can be obtained using the 
potential difference estimated above. However, even the 
extremely relativistic 100 MeV electrons which might be 
expected for a hot Jupiter will only penetrate down to 
about the 50 millibar level (Berger et at, 2005), which 
is clearly insufBcient. Some of this energy will conduct 
to lower parts of the atmosphere and interior, but it is 
unclear how efficient such processes might be, and thus 
energetic particle precipitation is unlikely to contribute 
substantially to interior heating although it will increase 
conductivity in the polar regions of the atmosphere. 

The second mechanism for energy transport I will con- 
sider involves the planetary global electric circuit (GEC). 
On Earth, the ground and ionosphere together form what 
is essentially a leaky capacitor with the atmosphere act- 
ing as a dielectric. The voltage difference of ^ 240 kV 
between the ground and ionosphere is maintained by a 
combination of thunderstorm action (« 80%) and the 
solar wind-driven potential (« 20%), and the circuit is 
closed through the fair-weather load resistance of the at- 
mosphere, which is approximately 250f2 (Tinsley, Burns, 
& Zhou, 2007). The total power involved in the ter- 
restrial GEC is thus of order 2.3 x 10^^ erg s~^, and 
can be an order of magnitude larger during substorms 
and other periods of enhanced solar activity. The ex- 
istence of a GEC on giant planets in the solar system 
is debated (for contrasting views, see Aplin 2006 and 
Simoes et al. 2012), primarily because the low conduc- 
tivity (cr < 10~" S m-\ Whitten et al. 2008) in the 
upper Jovian and Kronian atmospheres renders it prob- 
lematic; however, the higher temperatures anticipated in 
hot Jupiter atmospheres removes this obstacle. 

Although lightning may be available to generate poten- 
tial differences on hot Jupiters, as it is on Earth (Helling, 
Jardine, & Diver 2012), for simplicity's sake I consider 
only the FACs, which as noted above can be expected 
to produce a cross polar cap potential of as much as 
Ay = 320 MV between the day and night hemispheres 
on the canonical hot Jupiter with a semimajor axis 
a — 0.08 AU. This circuit can potentially be closed ei- 
ther via the ionosphere or the planetary atmosphere and 
highly conductive interior. However, the conductivity of 
the ionospheric route is dominated by the Pedersen con- 
ductivity of the ionospheric layer, which for reasonable 
estimate of ionospheric parameters (rig ~ 10^ — 10^ cm~'^, 
T ~ 10'^ — 10^ K) is roughly five orders of magnitude 
lower than the conductivity in the radial direction. Ac- 
cordingly, the resistance of the radial GEC circuit can 
be estimated using the calculated conductivity profile of 
Batygin, Stevenson, & Bodcnheimcr (2011, e = model 
as corrected on astro-ph) to be sa 1.7/iil. If we assume the 
canonical conversion efficiency of 1% seen in solar system 
magnetospheres (Vasyhunas & Song 2005), such a circuit 
leads to Joule heating levels of as much as 6 x 10^ erg s~^ , 
more than sufficient to sink the available magnetospheric 
power. Figure 4 shows the resulting heating profile for 
an evolved IMj, Tg = 1400 K planet in a form which 
can be directly compared to Figure 4 of Batygin, Steven- 
son, & Bodenheimer (2011); the profile from the GEC 
model is comparable to that of the Ohmic heating model 
at all pressure levels. Figure 4 also shows the bound- 
ary between the radiative and convective portions of the 
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Figure 4. A GEC heating model showing energy dissipation per 
unit mass, calculated for the conductivity profile corresponding to 
the evolved 1 Mj planetary model shown in Batygin et al. (2011). 
The black (solid) line corresponds to the radiative part of the plan- 
etary interior, while the red (dot-dash) line represents the con- 
vective inner region. Total heat input below the 1 kbar level is 
2.1 X 10^5 erg s"!. 



planetary interior, from which heating levels in the latter 
can be readily calculated. For the model shown, heating 
below the 1 kbar level is 2.1 x 10^^ erg s-\ or 0.06% of the 
intercepted stellar power Lp. This is entirely compatible 
with the required heating levels calculated by Burrows 
et al. (2007). 

Note that this model presumes a planetary Te = 
1400K; the effect of lowering this temperature to lOOOK 
is to raise the resistance by two orders of magnitude, 
which would lower the Joule heating rate by a similar 
amount, and which mimics the situation expected on so- 
lar system gas giants. The GEC heating mechanism is 
thus sensitive to the temperature structure of the atmo- 
sphere, but will generally be more efficient as semimajor 
axis a decreases (and thus atmospheric temperature and 
conductivity increase). One criticism of the proposed 
model is that the presence of an electrically insulating 
layer (as posited by Ohmic heating models) would sub- 
stantially reduce or eliminate GEC heating of the inte- 
rior. Such as layer might arise through the presence of an 
atmospheric inversion layer produced by opacity effects 
(Hubeny, Burrows, & Sudarsky 2003) or a thermoresis- 
tive instability (Menou 2012), and evidence for it at the 
10—100 mbar level has been detected in some, though not 
all, hot Jupiters (Knutson, Howard, & Isaacson 2010). 
However, it is likely that exoplanet atmosphere conduc- 
tivity profiles are spatially inhomogeneous; in particular, 
as noted above, nonthermal ionization due to auroral pre- 
cipitation will tend to increase conductivity in the polar 
regions that are most critical to the proposed model. 
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